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We consider two interacting Bose-Einstein condensates (EEC's) with different kind of the potential 
energy of interaction of the condensates: (a) the standard potential; (b) the potential has a positive 
three-body and a negative two-body scattering terms and (c) the potential has a positive four- 
body and a negative three-body scattering terms for the first BEC and a positive three-body and 
^— ■») , a negative two-body scattering terms for the second BEC. It is shown that in these cases there 

■ exist regular spherically symmetric solutions. Physically such solution is either a defect or a droplet 
created by the condensates. The defect is a cavity filled with one BEC on the background of another 
BEC. The droplet is an object on the background of the empty space. For (a) and (b) cases the 

^ . obtained objects are supported by a constant external trapping potential and for (c) case the droplet 

' is a self- maintaining object without any external potential. The possibility of construction of an 

' elementary logic qubit device on the basis of this droplet is discussed. 

(N 

^ ■ I. INTRODUCTION 

qh, Usually two interacting Bose-Einstein condensates are described by equations where the potential energy of inter- 
action of the condensates has terms of 2-th and A-th orders But in Ref. [H it was shown that the Hamiltonian 
\ of two interacting BEC's may have the interactions terms with both Q-th and A-th orders. There was shown that the 
^ 1 term of &-th order can be positive and the term of 4-t/i order can be negative. Here we would like consider three 
I ' different cases of two interacting BEC's: (a) both BEC's have standard interaction scattering terms and interact with 
■ O , a constant external potential; (b) both BEC's have positive three-body and negative two-body scattering terms in the 
^ ' potential and interact with a constant external potential; (c) one BEC is as in the previous item but another BEC 
has positive four-body and negative three-body scattering terms and they do not interact with any external potential. 

We will show that in such coupled equations set (describing two interacting BEC's) regular solutions appear. From 
our point of view such solutions describe either defect or a droplet created by two interacting BEC's. We will show 
that: in the case (a) we have a defect supported by a constant external potential; in (b) wc have a droplet trapped by 
a constant external potential and in (c) we have a droplet without any support of an external potential. The defect 
. means that we have a cavity on the background of the space filled with a first kind of BEC. The first kind of BEC has 
' a constant energy density at the infinity. The cavity is filled with another BEC whose energy density asymptotically is 
ps) . zero. The droplet means that we have an object filled with two BEC's with both asymptotically zero energy densities. 

■ In Ref. m it is shown that if three-body term is repulsive (the interactions terms of Q-th order in our language) and 
' two-body term is attractive (the interactions terms of A-th order in our language) then the conclusion is made: " Since 

psj ■ the two - body contribution to the ground state energy of a dilute Bose gas is negative, the three - body collisions 
in the regime where gs > could lead to the stabilization of the system. What is particularly interesting for such a 
system is that a boson droplet a boselet could become selfbound and the trapping potential is not required anymore 
to keep the particles together" . Below we will show that such self-maintaining configuration (droplet or boselet) does 
. exist as the solution of equations (p4|) ((35)) if one BEC has a positive four-body and negative three-body scattering 
^ ' terms and another BEC has a positive three-body and negative two-body scattering terms. 
- - - In this work we investigate the possibility of the existence of soliton - like solutions for two interacting Bose - Einstein 
condensates with and without any external trapping potential. The soliton- like solutions on the background of a 
trapping potential are widely discussed. In Ref. |3| the coupled Gross - Pitaevskii equations describing the dynamics of 
two hyperfine states of Bose - Einstein condensates are investigated and the integrability condition for the propagation 
of bright vector solitons is deduced. In Ref. Q the authors investigate the combined soliton solutions of two - 
component Bose - Einstein condensates with external potential. The results show that the intraspecics (interspecies) 
interaction strengths and the external trapped potential clearly affect the formation of darkdark, brightbright, and 
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darkbright soliton solutions in different regions. In Ref. Q the localized nonlinear matter waves in the two - 
component Bose - Einstein condensates with time- and space-modulated nonlinearities analytically and numerically is 
investigated. The dynamics of these matter waves, including the breathing solitons, quasibreathing solitons, resonant 
solitons, and moving solitons, is discussed. In Ref. sound waves in two - component Bose - Einstein condensates 
is investigated and proposed a new method of wave generation which is based on a fast change of the inter - species 
interaction constant. In Ref. Q stable skyrmions in two component Bose - Einstein condensates are considered. In 
Ref. [1] the authors report the numerical realization of robust 2 - component structures in 2d and 3d BEC's with 
non - trivial topological charge in one component. The vortex - bright solitary waves are found to be very robust in 
both in the homogeneous medium and in the presence of parabolic and periodic external confinement. In Ref. 0] a 
family of exact vector - soliton solutions for the coupled nonlinear Schrodinger equations with tunable interactions 
and harmonic potential is presented. 



II. DEFECT SOLUTION WITH TWO-BODY SCATTERING TERMS AND EXTERNAL TRAPPING 

CONSTANT POTENTIAL 

In this section we consider a defect solution. The solution describes the defect filled with one kind of BEC and 
placed in the space filled with another kind of BEC. 

Considering a two-component BEC, the behavior of the condensates that are prepared in two hyperfine states can 
be described at sufficiently low temperatures by the two-coupled GP equation of the following form [2[ 
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ih 



dt 
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where the condensate wave functions are normalized by particle numbers Ni — J 

)r the wave functions ipi.2{r,t) in dimension 

Y = (-V2+un|^/'i|' + ui2|^/'2|' + fi) V^i, 

(^-kV^ + U22 1-021^ + "21 IV'll^ + W2) V'2 



(1) 

(2) 

dV; Vi_2 are external trapping 



potentials. The resulting equations for the wave functions ipi.2{r,t) in dimensionless form can be written as 



dt 
.dip2 



(3) 
(4) 
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here k = mi/m2 and we redefined t/to — > t; r/Zo ~ — — , 

fl'<Pi (0) 



U12 



2mi(7i2^i''^0) 
W 



and Vj 



2miVi 



h^i^i (0); Vi,2 = j^] = ' 

For the simplicity we will consider the case A: = 1. We are searching for a 



static spherical symmetric solution: iJi.2{r,t) = tpi.2{r)- In this case Eq's ([3]) (j4]) are 

^1 + = ^1 [AaV^I + Ai (V? - M?)] ' 

+ ^ ^^^^^2 ^ (^2 „ ^2)] _ 



(5) 
(6) 



where Ai = ua^i = 1, 2; = — w;; A3 = U12 = "21- Eq's ([5]) ^ are written in the form which is convenient for the 
numerical calculations. For the numerical calculations we choose the boundary conditions and parameters values as 
follows 



Ai =0.1; 
^i(0) = l; 
^2(0) = VOJ; 



A2 = l.; 

^;(0)-0; 

02(0) = 0. 



X?. = 1. 



(7) 
(8) 
(9) 



The solution is searching as the nonlinear eigenvalue problem: /ii.2 arc eigenvalues and corresponding functions 0i.2 
are eigenfunctions. We solve equations set ([S]) ^ numerically. The profiles of functions V'1,2 hi Fig- [Dare presented. 
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The asymptotic behavior of the functions ■01.2 is following 



^1 



1p2 ~ -02 



(10) 

(11) 



where '01,2,00 are constants. 

One can consider Eq's ([5]) (El) as Euler - Lagrangian equations. Then the dimensionless energy density has the form 



eiM = ^(V0i)' + i(V02)' + V^(0i,2), 



(12) 
(13) 



Let us note that the energy functional is defined with accuracy of a constant. We choose the constant as — ^/^i — ^mI- 
In order to understand the physical sense of the obtained solution let us to introduce the energy densities for both 
EEC's 



2 , '^l /'_;.2 ,,2^ 
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Taking into account the asymptotic behavior (fTO|) ([TT|) we see that at the infinity 

Ai 4 



£(^-1,2) « -^Mi 



(14) 
(15) 

(16) 



i.e. is non-zero. It means that we have the space filled with BEC describing by ipi and at the center there is a defect 
filled with BEC describing by ip2- In Fig. [2] the profiles of the energy densities for both BEC's with tpi and ip2 are 
presented. 




FIG. 1: The profiles of ^1,2(2;), A*i = 1.61716, ^2 
1.49276 




FIG. 2: The profiles the dimensionless energy densities: 
curve 1 - the dimensionless energy density (|12p : curve 2 
- the dimensionless energy density of the BEC with -02, 
(|14|l : curve 3 - the dimensionless energy density of the 
BEC with ijjo. (fTSll. 



We see that choosing the external potentials Vi with some special values 



2mi 



4 



we have the solution describing a defect filled with one kind of BEC on the background of the space filled with 
another kind of BEC. From Fig. [2] we see the BEC filling the defect displaces the second kind of BEC. Formally it is 
analogously to Meissner effect in superconductivity. 

The profile for the dimensionless potential e(?/'i,2) as the function of functions ipi^2 hi FigI3]is presented. From this 
figure we see that it has two local and two global minima that gives rise to the regular solutions of Eq's ([Sj ([6]) . 




FIG. 3: The profile of dimensionless potential (|12[) as the function of ipi,2- One can see two local and two global minima. 



III. DROPLET SOLUTION WITH THREE AND TWO-BODY SCATTERING TERMS AND 
CONSTANT EXTERNAL TRAPPING POTENTIAL 



In this section we consider a droplet solution but BEC's have positive three-body and negative two-body scattering 
terms on the RHS of Eq's (fTS]) and (|T9| . 

In this case a two-component BEC can be described by the two-coupled equations (similar to CP) of the following 
form 
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dipi 
dt 



.,9-02 



2mi 
2m2 







4 






2 






2 




Uu 












KC/12 


^2 










4 






2 






2 




U22 


-02 




-V22 


-02 




hC/21 
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(18) 
(19) 



where the notations for Ni are the same as in Section |TT1 The resulting equations for the wave functions 01,2(^^7^) hi 
dimensionless form can be written as 



.901 

.di/j2 



( -V-- + Mil |-0i| -I'lllV'll +Wl2|V'2| -l-Wij'01, 
(^-kV^ + M22 1-02 I'' - W22 I V'2 1^ + M2I IV"! 1^ + '"^2) -02 



(20) 
(21) 



here k = mi/m2 and we redefined t/t^ t; f/lo to = ^^^L ; = "01 ^^^(0); V'1,2 = = ^"'^'^^t' — —\ 



2miVui,{'^{0) 



"12 = — — and «7i 2 



2mi Wi,2 



For the simplicity we will consider the case k — \. We 



are searching for a static spherical symmetric solution: Vi, 2(^,0 = i^i,2{r)- In this case Eq's ((20|) ([2T|) are 



^1 + -^i = ^1 [AsV-s + Ai {^l - M?) (3^'? - M?)] 



A30i' + 3A2V^2' ( ^1 - 3/^2 



(22) 
(23) 



where 3Ai = Uu; Vx = ifif; U12 = A3; 2A2^2 = ''^21; ^l^J■i = wi; W2 = 0- Eq's ([^ ([^ are written in the form which 
is convenient for the numerical calculations. For the numerical calculations we choose the boundary conditions and 
parameters values as ([7]) - ©. 

Again we search the solution as the nonlinear eigenvalue problem: /ii^2 are eigenvalues and corresponding functions 
■01^2 are eigenfunctions. The profiles of functions 0i_2 in Fig. 2] are presented. 



FIG. 4: The profiles of ^pi,2{x), ^ti = 1.0576294, ^2 = 
4.05682. 



FIG. 5: The profiles of the dimensionless energy densi- 
ties: the curve 1 - the dimensionless energy density (|26|) 
for both EEC's; the curve 2 - the dimensionless energy 
density ^ for ip2 BEG. 




FIG. 6: The profile of the dimensionless energy density 
(HHJ for ilJi BEG. 



The asymptotic behavior of the functions 'ipi_2 for such two interacting EEC's is following 



^1 



fJ-1 - V'l^oo- 



i'2 ~ 1p2. 

r 

where '01,2,00 are constants. 

One can consider Eq's ((22|) (p3|) as Euler - Lagrangian equations. Then the energy density has the form 



(24) 
(25) 



e(0i,2) = ^(VVi)' + ^(V02)' + y(0i,2) 



V (^m) = ^i^l {^l mlf + (y^l - ml) + '-^j^l 



(26) 
(27) 
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The energy densities for both EEC's are 



e2 = ^(V^2f + y^2MV^: 



(28) 
(29) 



Taking into account the asymptotic behavior (p4)) (|25|) we see that at the infinity 

e(V'i,2) ~ 0. (30) 
It means that we have a droplet filled with two interacting EEC's and trapped with the external potentials 

wi — Ai/i^ and W2 — 0. (31) 
In Fig's [5][6] the profiles of the energy densities for both EEC's and for ijji EEC and for ?/'2 EEC are presented. 

IV. DROPLET SOLUTION WITH FOUR AND THREE-BODY SCATTERING TERMS AND 
WITHOUT EXTERNAL TRAPPING POTENTIAL 



In this section we would like to consider the interaction between two EEC's where one EEC has a hypothesized 
positive strong four-body scattering term. The first EEC has three and two-body scattering (analogously to Section 
III[) but the second EEC has four and three-body scattering terms. In order to obtain a self -maintaining droplet 
without any external trapping potential we assume that there exists EEC with positive four-body scattering and 
negative three-body scattering. 

In this case EEC's equations are 



ih 
ih 



dt 

d'4'2 



2mi 

fi2 
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-C/ii 


^1 






+ C/l2 


■02 


1 



■01, 



U22 



■02 



V22 



05 



U2I 



01 



■02 



dt \ 2m2 

The resulting equations for the wave functions ipi,2ir,t) in dimensionless form can be written as 

.dijji 



' dt 
.d%l)2 



here k = mi/m2 and we redefined t/to 



(^-V^ + Wll |'0l|'^ - Ull 1-011^ + ■^11 1-011^ + U12 |V^2|^) 0'1, 
(^-/cV^ + U22 1-021^ - V22 |V'2|^ + "21 |'0l|^) ^02 

h - V-r'^'lO); 01,; 



(32) 
(33) 

(34) 
(35) 



2rnr<t^^2S^„^ ■ „ - 2miV-;°/^(0) ^, . _ 2mi-0y"(O) 
"11, Uii — 55 Uii, Vii — 7^2 



t; r/la f; <o 
Vii and Ui 



i/;i(0) ' 



Wll 



Vii aiiu uij = '^^^^j-^i — —Uij. For the simplicity we will consider 

the ease fc = 1. We are searching for a static spherical symmetric solution: V'i,2(f^, i) = V'i,2(?'). In this case Eq's ([34]) 
([55)1 are 



+ -Vi = V'l [AsV-i + Ai^? (V-? - t^l) (2^? - t^l)] 



V'2+-V'2 = 0-2 



(36) 
(37) 



where 2Ai = wn; 3Ai/.jf = un; Ai/i^J = un; A3 = U12 = 1*21; 3A2 = U22 and 2X2^1 = V22- Eq's ([M]) ([57]) are written in 
the form which is convenient for the numerical calculations. For the numerical calculations we choose the boundary 
conditions and parameters values as in ([T])-®. 

The solution is searching as the nonlinear eigenvalue problem: /ii,2 are eigenvalues and corresponding functions "01,2 
arc eigenfunctions. We solve equations set (|37p numerically. The profiles of functions ■0i,2 practically coincides 
with the profiles of the functions "01,2 from Section [IIII 
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The asymptotic behavior of the functions ■01.2 is following 



■01 ~ A^i - V'l.oo , (38) 

r 

02 ~ -02.00 (39) 

r 

where "01,2,00 SiTC constants. 

One can consider Eq's (j36|) (|37l) as Euler - Lagrangian equations. Then the dimensionless energy density has the 
form 

£(^1,2) = ^(V0i)Vi(VV2)' + F(0i,2), (40) 

V (^1,2) - ^i^t (0? - M?)' + Y^2 (^2' - f^l) + li^li'l (41) 



The energy densities for both EEC's are 

ei = ^(V^l)^ + ^VH^?-/^D^ (42) 



1 ' n2 , ^l.^,,4 /'„,,2 „2^ 



£2 = ^(V^'2)V^^|(^i-/.^). (43) 
Taking into account the asymptotic behavior ((38)) p9)) we see that at the infinity 

e(0i,2) « 0. (44) 

It means that we have obtained a self- maintaining droplet filled with two BEC. The profiles of the dimensionless 
energy densities pO)) . (|42p and (|43)) practically do not differ from the corresponding functions from Section Hill 

V. DISCUSSION AND CONCLUSIONS 

From the mathematical point of view regular solutions (having either a finite energy or at least finite values of the 
fields) for two interacting EEC's do exist only because the corresponding potential either p3)) or (l27l) or (|4T|) has local 
and global minima. If this potential has only a global minimum then following Derrick's theorem [l3| such solution 
may exist in the dimension l+I only I4| and this solution is known as a kink (Tsl - It is interesting to compare the 



properties of presented here solution with the kink solution. The main difference is that the solutions obtained here 
is topological trivial whereas the kink is topological non-trivial solution. The topological non-triviality means that 
there exists an integer (a topological charge) that does not change by a small deviation of the solution. Usually the 
topological charge is characterized by the solution behavior at the infinity. Let us emphasize once again that the 
solutions obtained here are topological trivial and they exist because the corresponding potential has both local and 
global minima only. 

Now we would like to discuss the physical interpretation of the obtained solutions. The interpretation strongly 
depends on the asymptotic value of the energy densities. In the case of non-zero value we have the space filled with 
one Eose - Einstein condensate 0i and on the background of this condensate there is a spherically symmetric defect 
(cavity) filled with another Eose - Einstein condensate 1^2- The asymptotic behavior of both condensates is different: 
for ^1 condensate the asymptotic value is non-zero but for ip2 condensate the asymptotic value is zero. The same is 
valid for both energy densities. If the asymptotic value is zero we have droplet solution: there exists a ball filled with 
two EEC's. There are two possibilities: either the droplet is trapped with an constant external potential (Section |III[) 
or not trapped without any external potential (Section lIVp . 

Numerical calculations show that most likely the existence of a regular solution depends on the potential form 
only: whether has it or not local and global minima simultaneously. In this connection one can mention that similar 
solution have been found for a scalar model of a glueball [l^- In this model the Lagrangian from SU(3) gauge theory 
by some manner can be approximated by a Lagrangian with two scalar fields. These two scalar fields approximately 
describe 2-th and A-th Green functions (i.e. the correlation between fields in two or four points in the spacetime) 
of SU(3) gauge fields. The first scalar field describes a gauge field belonging to subgroup SU{2) C SU{2>) and the 
second scalar field describes fields belonging to a coset SU{3)/SU{2). The kinetic terms from the initial Lagrangian 
gives rise to kinetic terms for both scalar fields. The terms like A^A'^A^A^ give rise to terms like {ipf — V'f-oo)^ 
in corresponding potentials for the scalar fields ipi (i = 1,2). Here Af^ is the SU(3) gauge potential. The difference 



8 



between the potential ([26]) and the potential obtained in Ref. [l^ is that the potential ((26)) is the polynomial of 
6-th order whereas the potential from [l^l is polynomial of i-th order. But in both cases corresponding potentials 
have local and global minima that leads to the existence of a regular solution. In Ref. a similar construction is 
considered and following conclusion is made "... a dilute Fermi - Dirac droplet will behave very much like a nuclear 
system." 

Finishing the comparison of obtained here solution with the solution obtained for the scalar model of a glueball we 
see that there is clear analogy between both solutions. It allows us to draw interesting and useful analogy between 
physical objects from high energy physics and EEC's physics: (a) the defect solution is similar to a cavity on the 
background of the space filled with a nonzero gluon condensate; (b) the droplet solution is similar to a glueball. 

Another important feature of the solution obtained here is following: In Ref. [Tsj it was shown that under certain 
conditions, the system of two EEC's has an almost degenerate ground state which is separated from the excited levels 
by an energy gap. The atomic ensemble then behaves like a two-level system that could be used to encode a qubit. 
The authors have offered the idea of using these two many-body states to encode a qubit and use it for quantum 
computation. One can use the defect /droplet obtained here as an elementary element (quantum gate) for creating 
the qubit device offered in Ref. [3l • 

Finally we would like to list the main features of the solutions obtained here: 

• The presented defect/droplet solutions exist for two interacting EEC's only. 

• The solutions exist for special choice of the EEC's potential: it should have local and global minimum only. 

• The solutions are eigcnfunctions for nonlinear differential equations. Consequently some physical parameters 
must have a well defined values (external potential, coefficients in front of n— body scattering terms and so on). 

• The droplet solution without any external trapping potential exists in the case of very strong interaction between 
atoms of EEC only. The droplet solution presented here does exist if there is a positive four-body scattering 
term of atoms for one EEC. 

• The defect solution exists for the standard potential of two interacting EEC's but with some special choice 
external trapping potentials. 
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